LOW LYING EIGENVALUES OF RANDOMLY CURVED QUANTUM 

WAVEGUIDES 
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Abstract. We consider the negative Dirichlet Laplacian on an infinite waveguide embedded 
in R 2 , and finite segments thereof. The waveguide is a perturbation of a periodic strip in 
terms of a sequence of independent identically distributed random variables which influence 
the curvature. We derive explicit lower bounds on the first eigenvalue of finite segments 
of the randomly curved waveguide in the small coupling (i.e. weak disorder) regime. This 
allows us to estimate the probability of low lying eigenvalues, a tool which is relevant in the 
context of Anderson localization for random Schrodinger operators. 
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1. Introduction 

In this paper we study eigenvalues of finite segments of randomly curved (quantum) waveg- 
uides. More precisely we derive lower bounds on the bottom of the spectrum of a waveguide 
segment of length L with perturbation of (a small) size k. The bounds are formulated in 
terms of the parameters L and k. In this context lower bounds are more challenging than 
upper ones. The waveguide consists of the set of points in M? which have a distance smaller 
than 7r/2 (half the width of the waveguide) to a random perturbation of a periodic curve in 
the plane. We consider the negative Laplace operator on this set with appropriate boundary 
conditions. As it has compact resolvent it has purely discrete, lower semi-bounded spectrum. 

Since the perturbation in each periodicity cell of the waveguide is determined by an in- 
dividual scalar random variable, the above bound can be considered as a multi-parameter 
variational problem. The main difficulty arises from the fact that the length of the waveguide 
is finite, but not fixed; indeed it tends to infinity. This also means that the number of varia- 
tional parameters becomes arbitrarily large. To compensate this, we have to chose the size k 
of the perturbations as a function of the length L. This is a special case of what is commonly 
called a weak disorder regime. 

The estimates on lowest eigenvalues we derive are motivated by the question whether a 
random Hamiltonian on a non-compact space has localized or non-localized states. One 
specific question in this context is: Is there a interval containing the bottom of the spectrum 
of the random Hamiltonian where the eigenvalues lie dense and where there is no continuous 
spectral component almost surely. The two first models where this question has been answered 
positively is the so called Russian-school-model (in [19]) and the Anderson model (in [1B S 
115]). In the later papers the so-called multiscale analysis was introduced as a method for 
proving (Anderson) localization, i.e. dense pure-point spectrum, with exponentially decaying 
eigenf unctions, almost surely. This method is an induction over increasing finite length scales. 
On each scale one considers a restricted Hamiltonian which lives on a finite volume. An 
important input for the multiscale analysis is the induction anchor, which is called initial scale 
estimate. It is this estimate which follows from our eigenvalue bounds, for the case that the 
random Hamiltonian is the negative Laplacian on a infinite randomly curved waveguide. The 
mentioned finite volume restrictions correspond to finite segments of the random waveguide. 

In a previous paper [3] we have studied the analogous problem for a different type of 
random waveguide, which we will call for definitness randomly shifted waveguide, cf. Section [J] 
for precise definitions and a comparison of the two models. Somewhat surprisingly, although 
the two types or random waveguides seem very similar, a very basic feature is different: 
namely, the configuration, which produces the lowest first eigenvalue. In one case the optimal 
configuration corresponds to zero randomness, in the other to maximal randomness. The 
model studied here is more difficult to analyse than the one of [3]. On the other hand it is 
also more natural, since it corresponds to the standard way how local curvature perturbations 
are introduced for waveguides, cf. e.g. [13] . 

An interesting feature of our model in comparison to other extensively studied random 
Schrodinger operators, e.g. the above-mentioned Anderson-model, is the geometric influence 
of the randomness. Another model with geometric randomness, which as recently attracted 
much attention, is the random displacement model, cf. e.g. [26[ [lj [18] [29] . One important 
feature of models with geometric disorder is that they exhibit no obvious monotonicity with 
respect to the individual random variables. This is the case for the random displacement 
model, for the randomly shifted waveguide (studied in [3]), as well as for the randomly curved 
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waveguide, to which the present paper is devoted. We refer the interested reader to the 
Introduction of j3], where we have discussed in some detail the challanges non-monotonicity 
poses in the spectral analysis of random Hamiltonians. 

There are also differences in the difficulties these models pose: for the random displacement 
model it is non-trivial to identify (all) optimizing configurations of the random parameters. 
(Here optimizing means: minimizing the first eigenvalue.) This is easier for the two mentioned 
models of random waveguides, although the optimal configurations turn out to be different 
in the two cases, cf. Section [3J On the other hand, in the random displacement model the 
random variables enter the potential, i.e. the zero order term, only, where in the waveguide 
models the higher order terms of the differential operator depend on the randomness. 

We would like to stress that localization for a certain type of random waveguides was 
established already in |25j. There the randomness enters via a variation of the width of the 
waveguide. This type of perturbation leads to a quadratic form which depends monotonously 
on the random variables and is thus structurally different from our model. 

As already mentioned, the main challange in our analysis is the non-monotone dependence 
of the Hamiltonian on the random variables. We have thus to overcome difficulties as they 
are encountered in other types of random Hamiltonians which exhibit a non-monotone depen- 
dence on the randomness. The model of this type to which most attention was devoted so far 
is the alloy-type random Schrodinger operators with single site potentials of changing sign, 
see e.g. |Z3>i3>0Z!>|2B!,0B!,|2D],|32]>0D!>|5D!- More recently also the discrete analog of 
this model was studied in [9], [32], [TO], [32], [33], [5], [8], [3JJ. Electromagnetic Schrodinger 
operators with random magnetic field [33], [33], [20], [31], [32], [3J, P3], p], [UJ, as well as 
Laplace-Beltrami operators with random metrics [37], [35], [36] exhibit a non-monotonous pa- 
rameter dependence which affects the higher order terms of the differential operator. Another 
relevant model (although not defined in terms of a countable family of random parameters) 
without obvious monotonicity is a random potential given by a Gaussian stochastic field with 
sign-changing covariance function, cf. [22], [45], [51j . 



2. Main deterministic results 

Let I > 1, I > a > and g £ Cq(R) ^ e an function with support in [0, a]. Let k ^ and 
p := (pj)j e z be a sequence of reals pj E [0, «]. We introduce the function 

G(-,p): R^R, G(x l3 p) :=J>;«&(a;i), 

j&Z 

where gj(x\) := g{x\ — jl) is the translate of the function g by jl to the right. Consider the 
curve Tp := {r p (t) \ t S R} C R 2 determined by the equation 

r p (t) := (t,G(t,p)), teM. 
By v p = v p (t) we denote the unit normal vector to T p at the point r p (t) defined by 

G'(t,p) 1_ 

P 1 1/2 (t,p) , P 1 1/2 (t, P ) / 
for tel. Here ' denotes the derivative w.r.t. t. We suppose that 

(2- 1 ) IMIc*[o,t] = L 

In the following when we say that a quantity does not depend on g or a, we mean that it can 
be chosen uniformly for all elements satisfying conditions (12. ip . 



V >® = I -^f?^— I ' - 1+ (G'(^)) 2 , 
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In a vicinity of T p we introduce new coordinates £ := (£1,^2) by the formula x = x(£) = 
r p(£i) + ^2^p(Ci)- In other words, given £1, we take a point at Y p corresponding to t = £1, 



Denote by T p (N,j) the upper and lower part of the boundary of D p (N,j), i.e., 

T p (NJ) := {x : jl < & < (j + iV)Z,£ 2 = 0} U {x : j7 < 6 < (j + = *■}. 

The remaining part of the boundary dD p (N, j) \ F p (N, j) is denoted by "f p (N, j). If j = 0, we 
will use the following shorthand notation: D p (N) := D p (N,0), T p (N) := T p (N,0), -y p (N) := 
7 p (iV, 0). Note that the length of D p (N) is L = N ■ I. 

In order to state the deterministic lower bound on the first eigenvalue, as a function of 
the parameters p = (pj)jez> we nrs t have to introduce a reference energy. This will be the 
spectral bottom of a comparison operator. More precisely, we will use the symbol H pcr (p) to 
denote the negative Laplacian in L2(D p (l)) with Dirichlet boundary condition on and 
periodic boundary condition on 7 P (1). Note that is this situation the array p = (pj)j consists 
just of the single value pq. So we identify here po with p. Let X pei (p) be the lowest eigenvalue 
of % per (/o) and ip per = ip peT (x,p) be the associated normalized eigenfunction. We extend the 
function ip peT ^-periodically w.r.t. £1 and denote the extension by the same symbol. 

We denote by H(p,N,j) the negative Laplacian on L2{D p (N,j)) with Dirichlet boundary 
conditions on T p (N,j) and with Robin boundary condition 



is the logarithmic derivative in ^-direction. We shall show in Sec. [7] that the function h 
is well-defined. In the following we denote for some s G M. by s ■ 1 the constant sequence 
Pj = s (j G Z). The important point is that ip per (-,p) is still an eigenfunction of the operator 
7~l(p ■ l,N,j) with the new boundary conditions, i.e. 



where we use for the restriction of ip per to the finite waveguide segment D p .j_(N,j) the same 
symbol. This relation was first exploited in the present context by Mezincescu [38]. For this 
reason we call (|2.2p in the sequel Mezincescu boundary conditions. We denote the lowest 
eigenvalue of H(p, N) = H(p, N, 0) by X(p, N). 

3 

Theorem 2.1. Assume a := ^ — a > 0. There exists a constant 5 = 5(a,g) > independent 
of p, I, N such that for 

(2.3) k < 5L~ n/2 , and < pi ^ k for all i G 0, ... N - 1 

the estimate 




(2.2) 




(0,x 2 , k) 



H{p-l,N,j 



KpW ev 



N-l 



(2.4) 




K 



3=0 
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holds true with 

J\\2 



r ,_ \}9 Wl 2 (o,i) ( 7T 16a , / 2 



^2 ^ y - ^2 \\y 11^(0,0 

In view of the well known result that curvature induces eigenvalues below the bottom of the 
essential spectrum of the straight quantum waveguide, cf.[II], it can be considered natural 
that maximal coupling constants produce the lowest first eigenvalue in the case of randomly 
curved waveguides 

3. Probability of low lying eigenvalues 

In this section we want to formulate an upper bound on the (small) probability that a 
randomly perturbed, finite waveguide segment has an eigenvalue close to the overal possible 
minimal spectral value. For this we use the deterministic results of perturbation theory, 
formulated in Section [2j The probabilistic results we obtain hold in the so-called weak disorder 
regime. To make this explicit, we use in this section a global (scalar) disorder parameter: For 
K>0we set 

ujj:=Pj/K (j G Z). 
Then we have the following equivalence for the condition used in Theorem 12.11 

Pj G [0, k] (j G Z) <==> Uj G [0, 1] (j G Z) 

We consider the parameters tOj,j G Z as a sequence of non-trivial independent, identically 
distributed random variables taking values in the unit interval [0, 1]. We denote the distri- 
bution of the single random variable tOj by pi, and by P the product probability measure 
on Q := [0, l] z which is the distribution of the random sequence tOj, j G Z. Of course O, is 
equipped with the natural product sigma-field. Now the deviations of the parameters from 
the extremal value, which determines the lower bound in Theorem 12. 1^ can be rewritten as 
follows: 

k — pj = k(1 — ujj) = k ■ Cj j for all j G Z 

Here we denote by ojj = (1 — ujj) G [0, 1] the flipped random variables, and the multiplication 
of a sequence to = (tOj)j^z with a scalar s G K by s • to. 

Now the result of Theorem 16.11 can be formulated in the new parameters as follows: If 

3 

a := — a > 0, then there exists a positive constant 5 independent of l,N,L and of the 
sequence to = {tOj)j^%, such that for all k < we have 



2 N-l 



(3.1) X(K-to,N)-X^(K) >C,6^-J>j whereC^^Ho'lH^a 

Here we denote the lowest eigenvalue of %{k-oj, N) = H(p, N) by X(k-lo, N) = X(p, N). Recall 
that A pcr (K) is the lowest eigenvalue of H per (p), i.e. the negative Laplacian in Z/2(-D p (l)) with 
Dirichlet boundary condition on r p (l) and periodic boundary condition on 7p(l). This is the 
lowest achieveable spectral value, if all perturbation parameters are bounded by k. 
The announced bound on the probability of low-lying eigenvalues is: 

Theorem 3.1. Let 7 > 22,7 £ N. Then there exists N\ = N±("f, pi,l,5, g,a) G (0,oo), such 
that for all N N± the interval 



I 



N 



C N~ 1/A ,5r 11/2 N~^ 
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is non-empty. For N ^ N\ and k G In, we have 

(3.2) pfwe n | X(k-u,N)- X pcr (K) < AT-I) < iV 1 "^ e-^LDP^/^ 

The constants C = C(fj,,l,g,a) and Cldp — Cldp(m) > 0; as we ^ as are given explicitly 
in Section^ 

Remark 3.2. The requirement k G Ijv encodes how our weak-disorder regime depends on the 
length scale N. For purpose of illustration, let us choose 7 = 33. Then 4 > 5^ = g, thus it 

is possible to choose k = const. N~±, which means that the allowed disorder regime does not 
shrink too fast for N — > 00. 

For the purposes of a localization proof, for instance using the multiscale analysis or the 
fractional moment method, the following initial length scale estimate or finite volume criterion 
is of interest. It can be derived from Theorem [3J] using a Combes-Thomas estimate [6]. This is 
a well established tool in the theory or random (Schrddinger) operators. An explicit derivation 
adapted to the setting of quantum waveguides can be found, e. g., in Section 4 in [3]. 

The mentioned initial length scale estimate is a bound on the probability that the Green's 
function exhibits exponential off-diagonal decay for energies in a small interval at the bottom 
of the spectrum, i.e. the overall minimum of the spectrum A per («). For the following statement 
note that for N large In C [0, 1]. 

Corollary 3.3. Let I, a, g, \l, In, 7, Ni and Cldp be as in Theorem \3.1\ Let Keljvfl [0, 1], 
a, P ^ 2 and set 

A := {x g D K . U {N) \ < xi < a}, 
B := {x G D R . U {N) \L-p^ Xl ^L}. 
Then there exists an absolute constant c such that for any N ^ N\ 

p(wgO|VA< X pci (k) + 1/(2VN): \\xa(H(k ■ u, N) - A) -1 xb|| < v 7 ^ e~ c ^( A > B )/^ 

4. Comparison with the randomly wiggled or shifted waveguide 

In this section we compare the results obtained in the present paper to those of [3], whether 
a related problem was studied. There we derived deterministic and probabilistic lower bounds 
on the lowest eigenvalue of finite segments of a randomly perturbed waveguide. In contrast to 
the model studied here, the random geometric perturbation was introduced by locally shifting 
or wiggling the waveguide. To explain the difference between the two models we recall the 
mathematical definition of the Hamiltonians studied in [3] . 

4.1. Statement of the results of [3] on randomly shifted waveguides. Similarly as in 
the model studied in this paper we consider random quantum waveguides in R 2 , determined 
by the following data: Let (u)k)keZ De a sequence of independent, identically distributed, 
non-trivial random variables with values in [0,1], k > a global coupling constant, I 1 the 
length of one (periodicity) cell of the waveguide, and g G Cq(0, I) a single bump function. 
The following function G:KxJl->K determines the shape of the waveguide 

G(xi,u) := ^uj k g(xi - kl). 
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Note that kG(xx,u) = G(x\, ku). For a global coupling constant k > , N G N and j G Z 
we define finite segments of a infinite waveguide 

D K)U (N,j) := {x G R 2 | j7 < x\ < (j + N)1,kG(x 1 ,u)) < x 2 < kG(x x ,u) + tt}. 

The upper and lower part of the boundary of D KjU1 (N,j), we denote by 

r K , w (iV,j) :={x G R 2 | j7 < X! < (j + N)l,x 2 = kG{x u lo)} 

U{x G R 2 | j7 < xi < (j + N)l,x 2 = kG(xi,u) +tt}. 

while dD KtUJ (N, j) \ T KtUJ (N,j) is denoted by 7« |(1 ;(iV,j). Denote the negative Laplace oper- 
ator on D K ^(N, j) with Dirichlet boundary conditions on T KyLJ (N,j) and Neumann b.c. on 
7K,uj(N,j) by H Kt w(N, j), and its lowest eigenvalue by X K ^(N,j). We use the following ab- 
breviation: 

D«,a,(iV) := 0), r«, w (JV) := T KjU (N, 0), 7«, W (JV) := 7*.u/(#, 0), 

H K ,UN) := 0), A«, W (JV) := A« |tl> (JV, 0). 

Note that D Ki(JJ (N,j) = Di tKW (N,j). Since k > is arbitrary we may assume without 
restricting the model 

(4- 1 ) maxdlflfUcp.q, ||||fl' / ||c[o ) q, ||ff"||c[o,q} = 1- 

Denote the distribution measure of lo^ by \i and by P = (^) fx the product measure on the 
configuration space Q = X/jg^fO, 1] whose elements we denote by u = (wfc)jfcgz- 

The minimum of the spectrum of the Laplacian on an straight waveguide segment Ho u{N,j) 
is equal to one, and no operator % Ki£J (iV, j) has spectrum below one. In this sense, one is the 
minimal spectral value for all random configurations. 

Theorem 4.1. Let 7 > 34. Let g and [i be as above and set 



1 f l U \A+ 3 lbllL 2 (o,z) 3 ll5 , ll| 2 (o,/) 
g{t)dt, c 2 = v , c 3 = ; 



9:9 I J 4 ' w ""' ~ z 2/ 3 ' ^ 5000 Z 7 

T/ien i/iere exists an initial scale Ni such that if N N\ the interval 

i_i 

2N-i 4 _ib 
,c 3 N z-f 



E{u k } 

is non-empty. If N N\ and k G In, then 

(4.2) P^weO X K ,U N ) ~ 1 < N ~^) < e~° hDvNlh 

/or i/ie constant Cldp > 0/ Lemma 15. .?), which depends only on fi. 

The corresponding deterministic lower bound on the first eigenvalue is: 
Theorem 4.2. Let I ^ 1, g : R ^ R, g £ R as above, and L := Nl. For 



8 



BORISOV AND VESELIC 



we have the bound 



N-l 



(4.4) 



X K AN) - 1 > 2 ll?lli a (o,0 



j=0 



4.2. Comparison of the results on randomly shifted and curved waveguides. For 

definitness we will call the model of [3| randomly shifted waveguide and the model studied in 
the present paper randomly curved waveguide. 

(i) In the first mentioned model the configuration of parameters which produces the min- 
imal first eigenvalue corresponds to a straight waveguide, where the perturbation is 
switched off. In contrast to this, the configuration of the randomly curved waveguide 
yielding the minimal first eigenvalue corresponds to a periodic waveguide with maximal 
possible curvature. This difference between the two considered problems was somewhat 
surprising for us, since, naively, the two models look very similar. 

(ii) Since for the randomly curved waveguide the bottom of the spectrum is achieved as the 
perturbation is maximal, the analysis of the model is more complicated. For we have 
to perform an asymptotic expansion not around the straight waveguide, but around the 
periodic waveguide %{p ■ 1, N) with maximal allowed curvature, which arises itself as a 
result of the perturbation. 

(iii) The smallness condition on the coupling constants — formulae (12. 3D and fj4.3f> — and 
the lower bounds on the first eigenvalue — formulae (|2,4p and (I4.4D — in Theorems 12.11 
and 14.21 exhibit a power-like, i.e. polynomial behaviour. The exponents do not coincide, 
but are of the same order of magnitude. 

(iv) In both models we employ a change of the variables x — > £ straightening the waveguide. 
After this change the operator depends on the variables pi. For the randomly curved 
waveguide the dependence on the p^s is not polynomial, but irrational. For the randomly 
shifted waveguide considered in [3] the relevant differential operator contains only terms 
which were constant, linear or quadratic functions of the variables p{. 

(v) In Theorem 12.11 we do not specify an explicit value of 5 as in the corresponding result 
in [3] on the randomly shifted waveguide. Although our technique allows us to obtain 
this constant explicitly, the required calculations are cumbersome and tedious. The 
main reason why the control of constants for the randomly curved waveguide is harder 
than for the randomly shifted one is the aforementioned irrational dependence of the 
perturbation on the random variables. 



We will use the following large deviations principle in the proof of Theorem 13.11 For a 
reference see for instance [7j. 

Lemma 5.1. Letuik,k G Z be an i.i.d. sequence of non-trivial, non-negative, bounded random 
variables. Then there exists a constant Cldp > depending only on p such that 



For the proof of Theorem{3j\we choose K G 2N and 7 G N and set N := IC . Thus N G N. 

1 1 

Set J = N/K = = N ~ . Following the same ideas as in |21|, 124"! [3] we decompose the 



5. Probabilistic estimates: Proof of Theorem 13.11 
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waveguide segment D K . UJ (N) into smaller parts 

(5.1) |J D K . u (K,j) 

i=o,...,j-x 

• 

where (J denotes a disjoint union, up to a set of measure zero. According to this decomposi- 
tion we introduce new Mezincescu boundary conditions on the interfaces joining the shorter 
segments (15. lj) . As in the original paper [38] we conclude that in the sense of quadratic forms 

J-i 

U{K-u,N)>Q)U{K-u,K,j). 

3=0 

If we denote the lowest eigenvalue of %{k ■ u, N, j) by A(k • u, N, j), it follows 

J-i 

(5.2) X(k-lo,N) ^ mm\(K-ui,K,j) 

3=0 

In particular, we have the inclusion 

J-i 

{wGfil A(k-w,A0 - A pc1 '(k) < Af~3 j C [J <jw G fi | \(n-u,K,j) - X pcr { K ) ^ K~l} 

3=0 

Since the random variables Uf^k G Z are independent and identically distributed, we obtain 
J-i 



| A(k- w,K,j) - X pcr { K ) < i^-i) TV 1 "^ p( w | A(k-w,K) - A pcr ( K ) < 

3=0 

For k < <5L~ n / 2 we have A(k • w, X) — A pcr (K) ^ Qftf^ S^o 1 ^3 an d thus the inclusion 

(2 



3=0 



1 V- 1 - <r ' 

Denote by E{u)o} = E{1 — wq} the expectation value of (any) Choose now n such that 



(5.3) -L* K-t ^ i.e. CX-i < k, where C ^ ' 



C^ 2 "2 ' ' yCi6E{w fc }" 

On the other we also need to satisfy the restriction k < 5(1 K)~ 11 / 2 = 5l~ 11 / 2 N^ 11 ^ 2 ^ . The 
upper and the lower bound for n can be reconciled if 7 > 22 and 



K ^ K x :-- 
The last inequality is equivalent to 

N > iVi 



4 

-22 



27 

2Z 12 



1 ■ \C lb E{u }5 2 
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For k satisfying ()5.3[) we are able to apply the large deviations principle of Lemma 15.11 and 
thus obtain 



1 K ~ 1 



3=0 



for K ^ K\. Hence, for N ^ iVi 

pjw | \(k-u,N) - X pci (k) < AT*} < AT 1 -^ e-^DpiW^ 

6. New parameters corresponding to a linearization around the optimal 

configuration 

Here we reformulate our main deterministic result, i.e. Theorem 12. II in a different parametri- 
sation, which reflects the perturbation theoretic proof. In what follows we want to study 
waveguide segments which correspond to configurations p where all coefficients pi are close to 
k. In other words, we perform an asymptotic analysis around a linearisation of the operator 
family, not around the point p = • 1 in parameter space, but rather around p = k ■ 1. Thus 
the natural parameters for perturbation theory are 

e := (ej) je z,tj ■= k - Pj 

In this and the following sections we denote by H K (e, N) the negative Laplacian on L2{D p (N)) 
with Dirichlet boundary condition on Y p {N) and with Robin boundary condition, cf. (12. 2p . 

d \ 
— -h(-,K)\ u = on i p (N), 

(6.1) where h{; n) : [0, vr] -)• R, 

1 dijj per 

is the logarithmic derivative in ^-direction. Recall that ■0 per ( - ? k) still satisfies the eigenvalue 
equation of the operator ~H K (0,N) with the new boundary conditions, i.e. 



?i K (0, N)ifjP ev = A(k)^ pct 

In the following we denote the lowest eigenvalue of H K (e, N) by A(e, k) and suppress here the 
^-dependence in this notation. Here is the announced reformulation of Theorem 16. 11 

3 

Theorem 6.1. Suppose a < There exists a constant 5 > independent of p, I, N such 
that for 

(6.2) k < 5L~ 11/2 , and p t < k for all i € 0, . . . N - 1 

the estimate 

/||2 / s N-l 



(6.3) A (e , «) - A-(«) > U-^)«E e ^ 



4L V 2 vr 



j=0 



ZioZds true. The value of 5 depends on g and a. 
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7. Deterministic lower bounds: Preliminaries 

In this section we present certain less technical preliminaries which are necessary for the 
setup of perturbation theory and the proof of Theorem 16. 1[ First we introduce the necessary 
notation. Given p, we define the operator 

T{p) : L 2 {D p {N)) -> L 2 (D (N)), (T(p)u)(0 := u(s(£)). 

corresponding to the change of the variables where £ is associated with p. If = k 

for all i, we will write shortly T(n) instead of T(k,k, . . .). Let us study in more detail the 
change x i->- £. It is easy to check that 



= AV f , 



A 



-l ._ 



/ dxi 


dx 2 \ / 




ik) = [ 








86/ \ 



■1/2, 



_ G /p-l/2 



-1/2 



We set 



^(£,p) :=1-6^1,P), Kfa,p) :=G"(^,p)P 1 3/2 (6, p) 
and obtain by direct calculation 

(7.1) P 3 := det^ 1 A = P~ 1/2 (l - £ 2 G"Pf 1/2 + G' 2 ) = P\ r ' l P 2 , 

where ' denotes the derivative w.r.t. £i- Hence, 

(7-2) II U Hl 2 (D p (7V)) = (T(p)u, PsT(p)u) L 2 (D (N))- 

Thus the operator 

7^r(p) : L 2 (D p (N)) -> L 2 (D (N)), (Vftr(p)«)(0 := V^(z(£)) 

is unitary. We observe that Pj(£, p) = 1 + O(k), i = 1,2,3, for p-,- < k. Thus for sufficiently 
small values of k, non of these functions vanish. One can calculate directly that 

(Pr 1 

P 3 A*A 



(7.3) 

(7.4) 

(7.5) 
(7.6) 





P 3 



|V r u 



*\\L 2 (D P (N)) =(AV € T(p)u, P 3 AV 5 T(p)n) L2(Do(7V)) 

2 



P, 



/2 dT(p)u 



+ 



L 2 (D (N)) 



P 



/ 2 dT(p)u 
96 



L 2 (D (JV)) 



A* 



1 



9 15 9 „ d 
+ —P3- 



\ =A i -Q p = n^(0) + Q p , 



p 3 \dZ1P3dh d& °d& 

d 2 „ d d 

26A-(6,p)-e 2 2 A' 2 (ei,p) 



G 2 (C,P) := 



1 9 1 



^(e,p) 



+ 



G'{£ x ,p)K£ uP ) 



1 <9P 3 _ Kfcup) 
Psdfr P 2 (tp)' 
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Under the standing assumption (|2.ip we can find a constant C independent of I, a, g, N, and 
k such that 

(7-7) || Sn \\ C (ptfN)) + II Ql Wc(D^(Nj) + II Q 2 \\ C (D^Nj) ^ Ck 

uniformly for all configurations satisfying pi ^ k for all i £ 0, ... N — 1. 

Our final aim is to have a lower bound on the movement of eigenvalues under a perturbation. 
However, first we need some rough a-priori information about the position of eigenvalues. 
This amounts of showing that the relevant eigenvalue is not too far from the one of a simple 
reference operator. So an upper bound on the distance is in question. For this preliminary 
considerations it will be sufficient to consider the operator on a waveguide segment of length 
I, i.e. the case N = 1. Before we state the next Lemma, we collect some simplifications which 
will be useful here and in subsequent considerations. 

Remark 7.1 (Simplified expressions for the case iV = 1). Here we consider the case N = 1, 
i.e. the operator T-L per (K.) on the unit cell D K (1). Note that in this situation there is only 
one coupling constant po in the game; for simplicity we will write here p = po and similarly 
<?(£i) = go(£,i). Furthermore the following functions take on a simple, explicit form: 

= i + (g"(6,p)) 2 = 1 + (pg'iCi)) 2 
k(Zi,p) = w"(6)^ 3/2 (6,p) = p</'(£i)(i + (pg'(ti))Y 3/2 
Pt&p) = i - = i - 6p/(£i)JT 3/2 (£i,p) 

p 3 (£,p) = Pi /2 p 2 = p1 /2 (Zup) - 6p/(6) J Pr 1 (6,p) 

o (f n\ K ^ p) elM 

In the following we will make use of the following error term estimates on the functions above: 

p} /2 (Zi,p) = 1 + y</(6) 2 + 0( P % ptHZi,p) = 1 - pV(6) 2 + 0( P % 

p 3 (e, p) = Pi /2 (6, p) - p&</'(Si)prH$i,p) 

= 1 - + y</(6) 2 + p 3 6</'(6)</ (6) 2 + 0( P A ). 

in particular Q p = 0{p). 

Lemma 7.2. There is a constant c independent of a, I, g such that for p < c the eigenpair 
X pcr (p), ip pcr (-,p) of the operator T~L peY (p) on the unit cell D p (l) satisfies the relations 

|A* er (p) - IK Cp\ \\T(pW ei (; P) ~ C " P^rWcHluT)) < °p 2 ' 

where the constant C is independent of p, a, I, and g. Here ipQ ei (£) '■= y/2/(irl) sin £2 "is the 
ground state o/'H pcr (0), and ip per is orthogonal to ipQ BT in L2(Dq{1)) and solves the equation 

(7.8) (H per (0) - = -fl^sp. 



SPECTRAL PROPERTIES OF RANDOMLY CURVED QUANTUM WAVEGUIDES 13 

Remark 7.3. The fact that the function ^ per depends only on the variable £2 leads to several 
simplifications. On the one hand, Q p ip peT = Q 2 (p) gf^cT 1 '■ ® n ^ ne other, integrals involving 
the functions g and ipQ can be often simplified by separation of variables, e.g.: 

W ^7— i^o er )L 2 (D (i)) = 7- / / d£ 2 cos(£ 2 )sin(£ 2 ) = 0. 

In the following this will be used repeatedly. 

Proof. The operator 7" (p)rl pct (p)T~ l (p) can be regarded as a perturbation of % per (0), since 
T(p)ri pcT (p)T(py 1 = H pcr (0) + Q p . Note that Q p is relatively bounded w.r.t. the operator 
"H per (0). This follows, since the boundary conditions of the functions in the domain of % per (0) 
imply that the boundary integral terms resulting from partial integration vanish. The bound 
(|7.7p implies that for sufficiently small p, Q p is relatively bounded w.r.t. rl pcr (0) with relative 
bound strictly smaller than one. Consequentyl the sum rl peT (0)-\-Q p is a closed operator on the 
domain of % pcr (0). For a ip in this domain and every 4> e L 2 (D (1)) the function p h-> (4>, Q p tp) 
is holomorphic for p in a small complex neighbourhood of zero. Thus p 1— )■ % pcr (0) + Q p is 
a holomorphic family of operators of type (A) in the terminology of [23} VII. §2]. Note that 
one is the lowest eigenvalue of % pcr (0) and ip per is the associated normalized eigenfunction. 
Thus, e.g., Theorem II. 5. 11 of [23] implies that there is a constant C such that 

A pcr ( P ) = 1 + P xr + p"xr(p), \xr( P )\ $ c, 

T(p)r eT (;P) = C r + + ✓VT M ll^ 2 per IIl 2 (d ()( i)) < C, 

where X peT 6 R and ip per : -Do(l) ~^ R are independent of /). The constant C is independent of 
p, a, /, and g since the perturbation Q p is uniformly bounded in these parameters, provided 
the normalization condition (12. It) holds. 

We substitute (17. 9h into the eigenvalue equation M per (p)ip peT (■ , p) = A per (p)V> per (-, p) and 
obtain 
(7.10) 

T(p){n pci (0) + Q p )(C r + PVT + pVH = H pei (pW ei (;p) = \ pcv (pW ci (;p) 

= up) (1 + P xr + P 2 Ar r ) « cr + p< cr + p 2 vr ) . 

Thus, 

H pc *m pcr + pu p ™mr + sx er + ps^r + ^(p 2 ) = c + p< er + pArvr + ^(p 2 )- 

We want to isolate the terms which are linear in the perturbation parameter p. For this reason 
we substract the eigenvalue equation % pcr (0)^Q er = ip per for the unperturbed operator, 

Bih per 

pw^mr + q p < cv + pQ^r = P n^m pci + q 2 & p)-^h+ o{ P 2 ) 

= p< cr + pA pc > per + O(p 2 ), 
divide by p, and finaly take the limit p — > resulting in: 



(^ pcr (o) - i)vr = -9"^- + Arvr 



This equation is solvable, if and only if the right hand side is orthogonal to ip p in L 2 (Do(l)), 
i.e., 

A pcr = Af r (v>r r ,< cr }L 2 (z)o(i)) = (g"-Q^-,i>o cr )L 2 (D (i)) = 



14 BORISOV AND VESELIC 

This implies (|7.8p . The orthogonality condition for ip pcT is implied by the identities 

i =ii^ per iii 2 (D P( i)) = (n P w ei ,p} /2 P2T(pw ei )L 2{Do{ i)) 

= (C r ^3C r )L 2 (Do(D) +2/>(C r ) ^r)L 2 (Do(l)) +0( P 2 ) 

=(^0"' C 1 - P^9i)^o CT )l 2 (Do(i)) + 2 p{^T ^T) l 2 {d {\)) + °(P 2 ) 

= l + 2p^^r)L 2{ D 0{ l))+O(p 2 ) 
since J^d^g"^) = 0. 

The functions T{p)i^ per and ip per belong to W2 1 (Dq(\)). By the standard smoothness 
improving theorems (see, for instance, [39^ Ch. IV, Sec. 2]) and by the smoothness of the 
function g we obtain that T(p)ip peT and ip per are also elements of W2 A (Dq(1)). In view of the 
embedding W 2 4 (D (1)) C C 2 {D^X)) [Ml Ch HI, Sec. 6] the functions T{p)^ pcT and ^> per are 
the classical solutions to the eigenvalue problem for tp per and to (17.8p . By applying Schauder 
estimates (see [Ml Ch. HI, Sec. 1-3]) we conclude that T{p)ip per and ip pcv belong to C 3 (Do(i))- 

It remains to prove that the asymptotics (|7.9p for T(p)tp per holds true also in C 3 (Dq(1))- 



norm. In order to do it, by (|7.10p we write first the equation for ip 

n pcr (o)^ per ={-Q P + A pcr )V2 er 



per 
2 



A per « cr + wn + A? c >r + p^Qp^r + P - 2 (q p - P9 "^A ^ 



+ T~\p) 

As above we again apply smoothness improving theorems from [39], [M] and obtain the 
uniform in p estimate 



per 




' P er I 



IIV2 IIC3(D () (1)) 

which yields the desired asymptotics for ip per . □ 

Since ip peT is the ground state, it is positive in D p {\). Hence, h in (|6.ip is well-defined. Let 
us check its behaviour near the Dirichlet boundaries. It follows from the definition of T{p) 
that 

{T( P w cr )(o,t2,p) = r cT (o,t2,p). 

In view of the smoothness of ip per we can apply Taylor formula as 6 — > +0, 

dT(p)ip peT 



^ per (0,e 2 ,p)=6- 



96 



2 9 2 T(p)V per 



6i=o 

«2=«2 



for some 6 6 [0,6]- Here we have also employed the Dirichlet condition for ip per at 6 = 0. 
We substitute the asymptotics for T{p)il> pcT from Lemma [7.21 into the last formula and arrive 
at 

V> per (xi,0,p) - \ — x 2 = 0(x 2 p + xl) as x 2 -> +0, p -> 0. 

V 7Tt 

In the same way we obtain 

< er (2;i,0,/))-yja; 2 = O((^-2;2)p+(^-X2) 2 ) as x 2 -> +0, p -)■ 0. 
Therefore, for small p this function satisfies the estimate 

(7-11) Maori < CP, 

where the constant C is independent of I, p, and 5. 
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After this preliminary analysis of the operator on the unit cell, i.e. a waveguide segment of 
length I, we turn now to waveguide segments composed of N cells. In what follows we denote 
all numerical constants by C. In different formulas this symbol will denote different positive 
numbers independent of L, N, I, p, a, and g. Let Hq(0,N) be the operator H K (0,N) taken 
for k = 0. In fact, Ho((i,N) is the Laplacian on Dq(N) with Dirichlet condition on Tq(N) 
and with Neumann condition on 70 (N). Note that the ground state of Ho(0,N) is a 
function of £2 only. Thus the function H defining the Mezincescu b.c. vanishes. 

The lowest eigenvalue of H K (0, N) is A per (fv); the corresponding eigenfunction is the function 
^ per . Let A per (K) be the second eigenvalue of H K (0,N). For the purposes of perturbation 
theory we need to establish a minimal distance between the two lowest eigenvalues A per (K) 
and A per (K) For this purpose we need 

Lemma 7.4. Under the assumption \2. 1}) there exists a constant C independent of l,N,a,g 
such that the eigenvalue X p obeys the inequality 

\\ p (k) - 1 - 4vr 2 L^ 2 | < Ck. 

Proof. The quadratic form associated with % K (0, N) reads as follows 



(M k (0,N)u,u) L2{Dk{n)) 



iVd 



L 2 (D K (N)) 



+ J h\u\ 2 dx2 — J h\u\ 2 dx2- 



xi=0 



x-i =L 



Denoting v := T(k)u, by (I7.4p . (I7.3P we rewrite the form as follows 



{n K (0,N)u,u) L2 (D K (N)) 



-1/2 



dv 



+ 



L 2 (D Q (N)) 



1/2 



dv 
W2 



L 2 (D Q (N)) 



+ I hfa,K)(H0,&W-\v(L,Z 2 )nd£2. 
lo 

= (Ho(0,N)v,v) L2 ( Do ( N )) + (Q K V,V) L2 (D (N)) 

Since the two lowest eigenvalues of ^o(0, N) are 1 and 1 + Air 2 L~ 2 , it is sufficient to estimate 
the quadratic form of Q K . It will be convenient to use the abbreviation (•, -)o := (•, •)l 2 (d (N)) 
and || • 1 1 := II " \\l 2 (D (N))i m the sequel. The estimate ()7. and standard embedding 
theorems yield 



(7.12) 

Thus we obtain 

(QKV,v) L2iDo{N)) 



h(t,K)(\v(o,t; 2 )\ 2 -\v(L,&\ 2 )d& 



< Ck|| V e v||o- 



+ J ^,«)(Ko,e 2 )| 2 -KL,6)| 2 ) 



d6 



+ 0(k\\Vv\\1) 
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=0 (k||Vv||§) 
Thus we have the estimate 

(H (0, N)v, v) + (Q K v, v) = (1 + 0(k)) (H o (0, N)v, v) 
Now we apply the minimax principle and obtain 

X p (k) = max min {(H o (0, N) + Q K )v, v) 

w£L 2 (D {N)) veW^(D (N)),v±w,\\v\\=l 

= (1 + 0(k)) max min (H o (0, N)v,v) 

w£L 2 (D (N)) ueW^-(D (JV)),w±it;,||t)||=l 

= (1 + 0(k)) (l + (^) 2 ) = 1 + (^) 2 + 0(«) 
since L ^ 1. □ 

The proven lemma implies that there exists C > such that the set H:={AeC:|A— 1|< 
Ck} contains no eigenvalues of H K (0,N) except A pc1 "(k). Now for A £ 5: 

\X p (k) - A| > |1 + 4^ 2 L" 2 - 1| + |1 + 4vr 2 L^ 2 - A p (k)| + |A - 1| 
^ 4vr 2 L" 2 - 2Ck 

Thus the distance from this set H to the remaining part of the spectrum a("H K (0, N))\{X P (k)} 
is estimated from below by 4-7r 2 L -2 — 2Ck. Theorem 16.11 concerns values of k which are much 
smaller than the inverse length 1/L of the waveguide segment. Thus we can certainly assume 
K < vr 2 /(CL 2 ) which means that 47r 2 L" 2 - 2Ck ^ 37r 2 L" 2 is positive. By [23l Ch. V, Sec. 
3.5] it implies that for A £ H the identity 

(7-13) (H K (0, N) - A) = >_ K ^ + K K (X) 

holds true, where tfj pcr := N~ 1 '" 2 ip per , (■,■}« := {',')l 2 (d k (N))i T^k is holomorphic w.r.t. to 
A € H as an operator from L 2 (D p (N)) to W 2 2 (D p (N)). The factor N~ 1 / 2 in the definition 

of ip per provides the normalization ||?/ ;Pcr ||L2(Do(A r )) = 1- For / E L,2(D K (N)) the function 
u = lZ K (X)f solves the equation 

(7.14) (n K (0,N) - X^( K ))u = I, f K :=f- ^(f, ft ei )L 2 (D K (iv)), 

(7.15) ||/ k ||l 2 (d k (7V)) < ||/||l 2 (d«(JV))- 

For each / G L2(D K (N)) the function TZ K (X)f is orthogonal to ipQ CT in L,2(D K (N)). 
The formulas ([72J> . (17131) imply 

(7.16) (TW^(0, iV)^ 1 ^) - A)- 1 = ^^In)-? ^^ 1 + ^ k(A) ' 

Let Ho(0,N) be the Laplacian on Dq(N) subject to Dirichlet boundary condition on To(N) 
and to Neumann one on 70 (N). 
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For the subsequent estimates in the perturbation theory we will need some a-priori bounds 
on several 'tame' operators, like resolvents. These are gathered in the next 

Lemma 7.5. For all A G S and all f G L2(D K (N)) the estimates 

(7.17) ll^(A)/||y^(D p (7V)) < CL 2 \\f\\ L2(DKm , 

(7.18) ||£«(A)/ - n (l)T(K)f\\ W2HDo{N)) < C7kL 4 ||/|| L2(Z)k (Ar)) 

ZioZd £rae ; where the constants C are independent of L, I, N , a, g, and f. 

Proof. Given and / G L,2(D K (N)), we introduce f K and u by (|7.14p . It follows from the 
definition of H and 1Z K (\) and (|7.15p that 

(7-19) \\u\\l 2 (d k (N)) < CL 2 \\f\\ L2{DK(N)) . 

Let x = x(t) be an infinitely differentiable cut-off function equalling one for t < 1/4 and zero 
for t > 1/2. We denote 

<H£, «) == x(£i)e^^ + X(L - tfeV-Wte'*') + 2 - x (&) - X (L - 
By (|7.1ip we have 

(7-20) U-l\\ c2{7m - }) ^CK. 

We construct u as u = T~ 1 (k)(/>u. We substitute this identity and (|7.5p into the equation 
(I7.14p which implies that for u 

(7.21) (H o (0, N) + C K )u = u + <t>- x T{K)f R , 

while C K is a second order differential operator such that 

( 7 - 22 ) \\£kV\\w 2 2 (D (N)) < C^IMIwa^-DoW)- 

It follows from (177191) that 

(7-23) II"IIl 2 (£>o(JV)) ^ C-^ 2 |l/llL 2 ( J Do(A r ))- 

Reproducing word by word the proof of Lemma 7.1 in [341 Ch. IE, Sec. 7] and employing 
(|7.23j) . one can prove easily an estimate 

(7-24) \\ v \\ W2HDo{N)) < C\\rl (0,N)v\\ L2{Do{N)) . 

Combining this estimate with (|7.22p . we can solve the equation for u as follows, 

u = Ho(P, N)-\l - £ K ^ (0, iV)" 1 )^ 1 ^ + ^T^/k), 

where all the operators are well-defined. Applying the estimates (|7.23p . (|7.24p once again, we 
arrive at the inequality 

H"llw2 2 (Do(A0) ^ CL H/lli2(£>o(iV)). 

which implies (|7.17p . 

We proceed to the proof of (|7.18p . We rewrite (|7.2ip as 

(7.25) (« (0, N) - 1)5 = <^- l T(K)f K + C K u. 

In accordance with (17.14p the function uq := Ho(l)T(n)f is the solution to the equation 

(7.26) (n (0,N)-l)u = f , 

fo ■■= T(«)/ - (T( K )/,< cr )L 2 (D oW )C r > C r == N- l ^o- 
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By the definition of 1Z K (\) and lZo(l), the function u = T 1 (k)4>u is orthogonal to ip pcr in 
L2(D K (N)), while uq is orthogonal to ipQ Cr in L2(Dq(N)). Hence, by the definition of (ft and 
Lemma 17,21 

(727 ^ u = u + ^ cr (u,^ cr ) L2{Do{N)) , 

{u^T)l 2 {d {N)) = 0, \(u^l cr ) L2{Do{N)) \ < Ck. 

The equations f?T26D imply 

(« (0, AT) - l)(n - n ) = <F X T{K)f K - C K u - /„, 

where (2 — uq) is orthogonal to ipQ er - Since the function u — uq is well-defined, the right hand 
side of the last equation is orthogonal to i^q CT in L2(Dq(N)). Moreover, the definition of (/), 
f K , /o, Lemma IT21 and (|7.22p yields 

||^ _1 T(«;)/k - C K U - foh 2 (D (N)) < C KL2 ||/llL 2 (D K (Ar))- 

Employing the aforementioned facts, by analogy with (|7.19p . (j7.24j) we get 

II s - u o\\w 2 2 (d (n)) ^ CkI? 11/11 l 2 {D k ( AT))- 
Since T(n)lZ K {\)f = T(k)u = u<j), the last inequality and (|7.27|) imply (|7.18|) . □ 



8. Deterministic lower bounds: Proof of Theorem 16.11 

After these preparatory lemmata we turn to the heart of the proof of of Theorem 16.11 It 
consists in deteremining the sign and estimating the three terms of the perturbation Q p = 

2n J|7 + Qi^j" + Q2^- In fact, since we have to perform a perturbative analysis not around 

the Laplacian "Ho(0, N) of the straight waveguide segment but rather around the one T-L K (0, N) 
with the maximal possible curvature, the effective perturbation will be "H K (e, N) — "H K (0, N) = 

Q P - Q K - 

As before we will work with the change of the variables i4( introduced above to transform 
the operator H K (e,N), 

(8.1) T{p)n K (e, N)T(p)' 1 = T(n)n K (0, N)T(k)- 1 +M(e, k) 

d 2 d d 

M(e, k) = M n — g + M 1 — + M 2 — , 
d£t 96 d& 

(8.2) M ll :=Q ll (^p)-Q 11 (^ K ), M t := p) - Q^, k) , z = 1,2. 

where the functions Qn, Qi have been introduced in (|7.6p . Let us calculate explicitly Q2 for 
values £1 G [0, 1]. First we expand the function p 1— > Qi{£,, p) around p = 0: 

(8.3) Q 2 (S,p) =g"^ 1 ) p (l + ^g"^ 1 )p + 0(p 2 )) = /(£l> + 6(/(£l))V + 

Suppressing for the moment the ^-dependence and using here the prime ' to denote derivatives 
w.r.t. the variable p we have 

Q 2 (p) = Q 2 (o) + Q 2 '(o)p + eV) 
Q , 2 , (p) = Q / 2(0) + O(p). 

Thus for some intermediate value p E [p, k] we obtain 

Q 2 (p) - Q 2 («) = Q! 2 ( K )( p - «) + Iq!'(k)( p - k) 2 + \q%{~p){p - «) 3 
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= -Q' 2 (0)e - Q£(0)eK + ^Q' 2 '(0) e 2 + 0( K 2 e) 



and after evaluating f|8.3[) and inserting the value at p = 0, 



follows, where we have made the dependence on £ explicit, again. Direct calculations of 
similar type show that 

(8.4) Mn := M§ } + m£ } + e^M^ , Mj := M^ 0) + Aff ) + e^ 2 Mf } , 

= 2 £j e 2 ^, m« = + Pj )(( 5 ;.) 2 - 3£ 2 2 (^) 2 ), 

^i (0) = ^2<, M« = - Ci ( K + Pj )(^ - 3£ 2 2 pf 

M 2 (0) = -e^J, M 2 (1) = - Ci (K + Pi ) 6(<#) 2 , 
for £1 G [(j - 1)1, jl], and 

(8.5) I M[f I + I M x (2) I + I M 2 (2) I C, 

where the constant C is independent of L, N, I, j, a, £, and g. 

Consider A £ 6. Employing the technique used in [2J Sect. 4], [T7] and (|7.13|) . (|8.ip . one 
can show easily that A = X(p) 6 is an eigenvalue of H K (e,N) if and only if it solves the 
equation 

A — X pcr (p) = {(l + M(e,K)li K {\))- l M{e,K)r c \P^), 

where ^ pcr := T(n)ijj pei , (•, •) := (•, -)o, P3 is taken for p = (k, . . . , k). We rewrite this equation 

as 

X(p) - X^(k) =(M(e,K)^ CI ,P 3 r ci ) - (M(e, K )K K (X(p))M(e, K )r cr ,P 3 ^) 

+ ((M(e, k)K k (X( P ))) 2 (I + M(e, k)^(A))^ 1 M(e, k)^ c \ P 3 r c "). 
The reduced resolvent 1Z K satisfies the resolvent identity 

- n K (fi) = (A - p)n K (p)1z K (x). 

Substituting this identity with A = A(p), p = 1 and (|8.4p into (|8.6p . and using the fact that 
72 K (1) = 72.o(l)) w e obtain 

(8.7) X(p) - A pcr («) = 5i(e, k) + (A(p) - A pcr (K))S 2 (e, k) + 5 3 (e, «), 

Si(e,«) = (M(e, K )^,P 3 ^) - (M(e, K )K (l)M(e,K)r ci ,r ci ), 
S 2 (e,K) = -(^(e,K)7e K (AP cr ( K ))7e K (A(p))^(e, K )^ er ,P 3 ^ per ) 
S 3 (e,K) = -(M^He, K ){l-n K (X^(K))M(e,K))r c \P 3 r ci ) 

- (M(e, k)TZ k {X^{k))M^ (e, K )r cr , p 3 rn 

+ ((M(e, K)H K (X(p)))\l + M(e, k)tL(A))- 1 M(e, k)^,P 3 ^), 

- (M(e,K)(U K (X(p))-n (l))M(e,K)r er ,P 3 rn 

+ (M(e, K )r c \(P 3 -P 3 )r cr ) 

- (M(e,K)K (l)M(e, K )r c \ (P 3 - l)^ er ), 
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.M(e,/c) :=7W (0) (e,K) + 7W (1) (e,K), 

M®(e, K ) := + Mf } A + M 2 W * 4 = , 1, 2. 

<9£f <9£i <9£ 2 

By (|8.5j) . (|6.2p and Lemmas I7.2( 17.51 we get the estimate for ^(e, k): 

/ jv \ 1 / 2 

(8.8) |5 3 (e,K)| < CK 2 \e\ 2 L 9 / 2 , where |e| 2 := ( Yj^J 

It follows from (|2.ip . (|6.2p . the definition of Ai(e, k) and Lemma 17.51 that 

|5 2 | Ck 2 L 2 sC CcL- 7 . 
Hence, for S small enough we have IS2I < 1/2, and by (|8.7|) . (|8.6|) we get 

(8.9) X(p) - X^( K ) > ^£^1 _ CK 2 \e\ 2 L 9 ' 2 . 

In order to control the contribution S\ we split it into two parts, using thereby identity ([8 
and Lemmas 17.21 

(8.10) 5i(e,«)=54(e,«) + 5 5 (e,«), 

s 4 (e,K) -jyr 1 ^, k^.O + ^ 1 (-M (0) < cr ,C r ) 

+ ^-^M^^r^D - iv- i (^(°)(6,K)c r ,6G'Vo pcr ) 

-iV- 1 (^( )(6,^o(l)^ (0) (6, K )C r ) C r >> 

S 5 (e, «) :=Si(e, k) - 5 4 (e, re), 
where G" is taken for p = (re, . . . , re). By Lemma 17.51 and (]8.5[) we can estimate £5, 

(8.11) \S 5 (e,K)\ ^CK 2 \e\ 2 L 7 ' 2 . 

We shall show that S5 can be regarded as a small order perturbation of the main term 54. 
In order to it, we need to calculate S4. It is easy to check that 

Ho"!! 2 N 

N-\M&k)1%*,1$*) = 2L E ej ' (K + Pi) ' 

8.12 

llo"ll 2 N 
-N-\M^{e,K)i>r^G"^) = - '^ 3(0,0 «X)6 J -. 

i=i 

By ([831) we obtain 

(^(°)( e , K )< cr ,c r ) + (^ (o) (e,^c r ^r) 

d 2 ih pei 8ib per 8ib p( 



/oC , c W^T /Per 



j=l \ ^ 2 / L 2 (D (j-l,j)) 
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d 8 d & ^2 /L 2 (Do(l)) 



i 2 (Do(l)) / j=i 

We take into the account the equation (17, 8j) for ^j 61 " and integrate by parts, 
(8-13) 

//^fyf* t ///^r 1 c //^r" ,per\ /.per c „ " \ 



, n2 / per / A \ 

/ /- // u rl / per \ t- II i per / per \ 



1 1 9" 1 1 L 2 (0,0 / „ .per 



11,2(0,0 

2/ ~V ri ' <9£ 2 /ia(A)(l)) 



Now consider the last contribution to S4. Denote u := 7£o(l)-M(°)(e, k)^o- This function 
solves the equation 

JL r)?/; per 

(8.14) (H O (0,N) - l)u = -g- 

and is orthogonal to tpQ er in L2{Dq{L)). We use this equation and proceed as in (I8.13D . to 
calculate 



per\ 



^ \ 3 <9£i J <9£i J <9£ 2 /L 2 (DoO-i,j)) 

/* c ;Per\ 

E^(^'(-^'^ ! - + S + #+«).<") , , „ 

ll«"ll 2 ^ W fl / Per 

119 IIl2(0 '°Es 2 + 2 E^> 



2/ f-f 3 f-' ' ^2 /i 2 (A0-U)) 

iV ^ 



22 BORISOV AND VESELIC 

We substitute the last identity and <KT2\i . <KT3\i into (pTL0l) . 



96 ' L 3 (D (j-l,j)) 
3 =1 

Note that terms involving ||</'||l 2 (o,;) cancel. Next we study each term in the obtained ex- 
pression for 1S4 . 

Lemma 8.1. The inequality 

fnir\ at \ ^ H^lliaCO,!) ( 7T 16a \ 

(8.16) 5 4 (e,/e)> I - - — r J «2^ej 

ZioZds true. 
Proof. Since 

(8-17) -7— (6,7T-6) = -^T- 6,6), 

the solutions ^f™, ii to the equations (|7.8p . ()8.14|) satisfy the same identity, 

(8.18) vr(6,vr-6) = -< er (6,6), = -^(6,6). 

It implies 

AT ^ ^ N Q 

^ ^ per 



(8.19) 



96 ' i*t(p Q (j-i,j)) 



o , per A^ 

2iv- i «/^ er ,/^_\ y, 

V 1 ' y 8Eo /Lo(D n (l))^ 



96 /l 2 (a,(i))^ j 

Let us estimate the first term in f|8. 18|) . We integrate by parts taking into consideration 
(EUD, (HH]), 

fri ej \ u,9j 96 /i 2 po(i-i,i))"^ e Aser*' 96 /wDoQ-ij)) 

3=1 3=1 
Vl|2 AT o N 

3=1 ^ 3=1 
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We construct the function u by the separation of variables, 

d^T _ ^ AT 4m 



a m sin2mg2 ; ol 



Q ~m Vvr/4m 2 -l" 

m=l 

oo 

(8.21) u(0 = E a m [/ m (&)sin2m£i, 

m=l 

where the functions U m solve the boundary value problems, 

N 

-U^ + (Am 2 -l)U m = -Y,^9j m L 2 (0,L), U' m (0) = U' m {L) = 0. 

i=i 

It follows from (|QJlt . (fg^TD that 

AT „ ,per Ho'll 2 AT oo AT 

(8.22) = + «£ •»«-(£«>*.«.) 

j=l ^ j=l m=l j=l 

We represent J7 m as 

N 

Um=7, £j9j + (4m 2 - l)U m , 

3=1 

where U m is the solution to the problem 

AT 

-U^ + (4m 2 -l)U m = -Y,^9j m L 2 (0,L), ^(0) = U' m {L) = 0. 

i=i 

Thus, 

at 

II ^» Ilia (0,0 + ( 4m ' " !) II 1113(0,1) = -(E e J*''^ 
AT AT AT 

( E Wj,U m ) L2{ol) = \\g\\l m E S 2 + ( 4 - 2 - 1) ( E e i*i> % 
j=i j=i j=i 

N N 

(4m 2 - l)||C/ m |U 2(0 ,L) < E e J* r , n n = 11^11^2(0,0 E e i- 

"i=i L2( °' L) j=i 

Two last relations imply 

N 

■ ^ /L 2 (0,L) 
and by (|g^2l it yields 

(8.23) f>(^#) ^Wkijfj 

3=1 J=l 

We proceed to the second term in the right hand side of (|8.18p . We construct ipf ev as 

56 V ttZ 



La(0,L) 
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where ^ cr solves the boundary value problem 

(8.25) (A + l)^ er = in £> (1), ffi* = g as & = tt, = -g as 6 = 0, 

and satisfies periodic boundary condition on £i = and £i = I. Hence, 

tf-ytofil) = 11^111(0.0 _ JX/ff* ti f*!F 



j 2g N <9& / £a(A>(l)) i V 7rZ \ 96 / La(A)(l)) 



l \ irl\ d£i d& ' l 2 {d {\)) 

Now we construct ^ er by the separation of variables, 

m=l ^ ' 

where (3m^ are introduced as the coefficients of the Fourier series for g', 

(8.28) g'(^) = f; cos + *4 _) ^ • 

m=l ^ ' 
The functions are defined as the solutions to the boundary value problems 

— l\V m = in (0,7r), y m (0) = -l, y m (vr) = l. 

The functions can be found explicitly, 



sinhA m (6-f) . /4vr 2 m 2 " 4vr 2 m 2 



sin 14 

v.(60-^.-i, if ^-1. 

We substitute (|g37j) and (18^8]) into (|5^B|> . 

,, 31) _ _ >mm + 1 (WL+)) . + (/4 -, )>n 



m=l 



where 



4^ m coth ^ ., 47r 2 m 2 
m " -(A^ + l) ' * Z 2 >X ' 
(8.32) B ._ ^ m cot^ 4vr 2 m 2 

r -— L if 4?r2m2 - 1 

£>m • 9) 11 ,9 

7T Z Z z 

In view of the estimates 

4zcoth^ 4 4zcot^ 8 , , 

' z^O, _ 1= < - T , 2 6 [0,1, 



7T 



V^TT tt' vr\/r^ 2 " vr 2 
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the coefficients B m can be estimated as follows, 

21 „ I Al . I 
> B m > 5 — , if 2m > -, > B m > 5—, if 2m < -. 

We substitute these estimate in (|8.3ip . 

_ 9 / 7/)P er ,,^ pe \ > %ik M1 _ 4/ v (/j+y+^-y 

V 1 ' 5 56 /L 2 (A,(l)r / VT3 2^ 

2m^-i- 

7T 2 ' m 

2m>l- 



2 llg / Hi 2 (o,0 _ 2Z ^ (A + T + (^T 

Z 7T 2 ' 771 



m=l 

Together with (|8.23[) . (|8.19p and an obvious estimate e,- ^ k it yields 

.33) 5 4 (e,«) > [ || L2(0i/) " ^ 2^ ZZ^ 

m=l / j=l 



It remains to estimate the last sum in this expression. 
We introduce the function 



m=l ^ ' 
I 

K--=\ J 3(6) d£i, 


and by the definition (|8.28[) of the coefficients f3m^ we see that 

wiLm = 5 £ (i) (WL +> )' + (A: 1 ) 2 ) + (4 +> ) 2 '. 

771=1 ^ ' 



2vr/3 Q 



^ 4 J ll«llLa(0,I) - 2 Z> + j J - -72-11^11^2(0,0 



m=l 



i.35) — > = — )l 2 ( ,z). 

m=l 



Now we use the fact that suppg lies in a segment of the size a and an obvious estimate 

7T 
ft 



ll<7 IU 2 (o,/) ^ — H<?IU 2 (o,z)> 

to obtain 



Z> > -^2ll«llLa(0,I)ll» llLa(0,I) 

m=l 



16 1, n n /u 16a.. ,|| 9 

> -— ll#IU 2 (o,z)lls IIl 2 (o,/) ^ h 2 (o,i)- 
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We substitute this inequality into (|8.33|) and complete the proof. □ 
We substitute the inequalities f|8 . 1 6 j) . (|8. 10j) . (jHTTTj), and a trivial one 

N 

\e\ 2 ^J2 e i 

3=1 

into (pOj) . 

aw - r W > (^|^ (I " ^) " CKi9/2 ) *X> 

Choosing the constant 5 in (16.2)) small enough, we arrive at the estimate (16. 3|) . 
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